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STATEMENT  OF  THE  PROBLEM 


Recently,  a  new  type  of  spatial  soliton\  based  on  the  photorefractive  effect,  has 
been  predicted  and  observed  in  both  a  quasi  steady-state  regime^’^  and  more  recently  in 
the  steady  state  regime."^^®  Compared  to  Kerr  spatial  solitons,  the  most  distinctive 
features  of  photorefractive  spatial  solitons  are  that  they  are  observed  at  low  light 
intensities  (in  the  mW/cm  range)  and  robust  trapping  occurs  in  both  transverse 
dimensions.  Both  of  these  attributes  make  photorefractive  solitons  attractive  for 
applications  and  for  fundamental  studies  involving  the  interaction  between  spatial 
solitons.  ‘  In  this  project  we  have  investigated  the  formation  of  nanosecond  high- 
intensity  self-induced  waveguides  in  order  to  produce  materials  with  enhanced  parametric 
conversion. 

One  transverse  dimension  theory  of  photorefractive  screening  solitons^’’  predicts 
a  universal  relationship  between  the  width  of  the  soliton,  the  applied  electric  field,  and 
the  ratio  of  the  soliton  intensity  to  the  sum  of  the  equivalent  dark  irradiance  and  a 
uniform  background  intensity.  We  refer  to  this  curve  as  the  soliton  existence  curve. 
This  existence  curve  is  important  because  experiments  show  that  considerable  deviations 
(~20%  or  more)  off  the  curve  lead  to  instability  and  breakup  of  the  soliton  beam^°’^^, 
while  much  smaller  deviations  are  typically  tolerated  and  are  “arrested”  by  the  soliton 
stability  properties.  In  the  case  of  a  low-intensity  photorefractive  soliton  beam,  i.e., 
intensities  in  the  mW/cm  to  KW/cm  range,  recent  one  dimensional  experiments  have 
shown  good  agreement  with  this  universal  relationship.^'’’^” 

Although  the  low-intensity  feature  of  photorefractive  spatial  solitons  is  attractive 
for  applications,  high-intensity  (MW/cm^  to  GW/cm^)  photorefractive  solitons  are  also 
interesting  since  the  speed  with  which  the  steady-state  screening  soliton  forms  is 
inversely  proportional  to  the  optical  intensity.  As  we  show  below,  solitons  in  SBN 
(strontium  barium  niobate)  can  form  at  nanosecond  speeds  for  GW/cm^  intensities.  This 
implies  that  for  photorefractive  semiconductors,  which  have  mobilities  100-1000  times 
larger  than  that  of  the  photorefractive  oxides,  soliton  formation  should  occur  at 
picosecond  time  scales  for  similar  intensities.  For  these  intensities,  however,  the  excited 
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free  carrier  density  is  no  longer  smaller  than  that  of  the  acceptors  and  the  space  charge 
field  is  due  both  to  the  free  carrier  and  the  ionized  donor  contributions/ 


In  this  project  we  report  the  first  experimental  observation  of  high-intensity 
screening  solitons,  along  with  a  comparison  between  experimental  results  and  theoretical 
predictions.  To  be  in  the  high  intensity  regime,  one  must  satisfy  the  requirement  that  1/r 


«a(uo^+l)«l, 


(Idark  fc)  . . 

where  r  -  Na/NA ,  a  -  s  r-  ,  Na  the  total  donor  number  density, 


Na,  the  number  density  of  negatively  charged  acceptors  that  compensate  for  the  ionized 
donors,  Uo^  is  the  ratio  of  the  soliton  intensity  to  the  sum  of  the  dark  and  background 
intensities,  Ib  the  background  intensity  (used  to  control  the  effective  dark  carrier  density), 
Idark  the  dark  intensity,  s  the  photoionization  cross  section,  and  y  the  recombination  rate 
coefficient.  In  our  case,  bright  high  intensity  solitons  in  photorefractive  SBN  can  be 
realized  at  incident  intensities  of  the  order  of  100  MW/cm^  on  a  background  of  Ib  ~ 
lOMW/cm^,  which  results  in  free-electron  densities  (for  T~300K)  of  10*^  cm'^  in  the 
center  of  the  soliton  and  10^^  cm'^  far  from  the  center.  For  the  experiment  reported  here, 
we  have  used  a  crystal  with  Na ~  lO^^cm'^  and  Na~  lO^^cm'^  (r=10^  and  a  ~  0. 1). 


For  the  experiment  a  Q-switched  YAG  laser  is  used  to  generate  a  high  intensity  8  ns 
second  harmonic  pulse  at  530-nm  which  is  split  into  two  beams.  One  beam  acts  as  the 
soliton  beam,  while  the  second  beam  provides  the  background  intensity.  The  background 
beam  fills  the  entire  crystal,  while  the  soliton  beam  is  focused  with  a  cylindrical  lens  to 
11  |im  (FWHM)  X  2mm  at  the  crystal  entrance  face.  The  cylindrical  lens  is  used  to 
observe  a  one-dimensional  (1-D)  soliton.  That  is,  the  beam  is  essentially  infinite  in 
extent  in  the  vertical  or  b-direction  and  1  Ip  in  the  horizontal  or  c-direction.  Both  beams 
propagate  along  the  crystalline  a  axis  (for  SBN  a  =  b)  and  the  electric  field  is  applied 
along  the  c  axis.  The  soliton  beam  is  extraordinarily  polarized  in  order  to  take  advantage 
of  the  crystal  large  r33  electro-optic  coefficient,  while  the  background  beam  is  ordinarily 
polarized,  making  it  possible  to  select  only  the  soliton  beam  for  observation  by  using  a 
polarizer.  In  particular,  the  soliton  beam  is  observed  using  an  imaging  lens  to  image  and 
magnify  the  intensity  distribution  at  the  crystal  entrance  and  exit  faces.  Fig.  1  shows 
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typical  experimental  results  of  the  entrance  intensity  waveform,  the  exit  intensity 
waveform  without  trapping  (zero  voltage),  and  the  exit  intensity  waveform  with  trapping 
(when  1500  volts  are  applied  to  the  crystal  between  electrodes  separated  by  6mm).  The 
waveforms  shown  evolved  using  8  ns  pulses  at  an  intensity  of  lOOMW/cm^  after  10 
pulses. 

SUMMARY  OF  THE  MOST  IMPORTANT  RESULTS 

While  it  takes  10  pulses  to  reach  the  steady-state  condition,  the  intensity  of 
illumination  during  each  pulse  is  in  the  high  intensity  regime  and  the  beam  diameter 
reaches  e'^  of  its  steady-state  value  during  the  first  pulse.  Moreover,  the  crystal  dark 
current  is  so  small  that  the  measured  beam  diameter  at  the  end  of  each  pulse,  was  found 
to  be  independent  of  the  time  between  pulses  over  a  range  from  0.1s  to  1000s.  For  these 
two  reasons,  our  crystal  can  be  considered  to  be  in  effective  illumination  at  100  MW 
intensities  for  a  period  of  80  ns  before  steady-state  is  reached. 

One  dimensional  high  intensity  screening  solitons  obey  the  normalized  nonlinear 
wave  equation^ 


d^u/d^^  +  5u  -  u/(l+u^)^^  =  0 


(1) 


where  u(^)  is  the  soliton  amplitude  (as  a  function  of  the  transverse  coordinate  ^)  divided 
by  the  square  root  of  the  sum  of  the  background  and  dark  intensities,  Uo  is  u(^=0), 
5  “2[(uo+l)  -l]/uo^  and  4=x/d  where  we  have  the  following:  d=(k^nb^reffV/^)'^^,  k 

=2nb/A,,  X,  is  the  free-space  wavelength,  nb  is  the  unperturbed  refractive  index,  reif  is  the 
effective  electro-optic  coefficient  for  the  geometry  of  propagation,  V  is  the  applied 
voltage,  and  ^  is  the  width  of  the  crystal  between  the  electrodes.  Equation  (1)  can  be 
integrated  numerically  to  obtain  the  spatial  profile  of  the  soliton  and  the  full  width  at 
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half-maximum  (FWHM)  of  the  intensity  as  a  function  of  uq,  which  is  the  soliton 
existence  curve.  These  results  are  shown  in  Fig.  2  (solid  curve)  for  the  range 
O.l^uo^lOO,  along  with  the  low  intensity  case,  (dotted  curve),  for  comparison.  The 
difference  in  the  two  theoretical  curves  is  because  the  change  in  the  refractive  index.  An, 
is  proportional  to  (1+u^)’^^^  for  the  high-intensity  solitons  of  Fig.  2,  whereas  for  low 
intensity  screening  solitons  An  is  proportional  to  (l+u^)’\  Our  experimental  results  with 
ID  solitons  of  the  same  width  and  wavelength  and  in  the  very  same  crystal  for  both  low 
and  high  intensity  solitons  are  marked  by  the  open  and  filled  squares,  respectively,  in 
Fig.  2.  It  is  apparent  that  there  is  good  agreement  between  experiments  and  theory  for 
both  the  high  and  low  intensity  cases.  The  predominant  reason  for  the  discrepancy  is  that 
the  background  beam  is  slightly  guided  by  the  refractive  index  change  induced  by  the 
soliton**^  (since  ri3  is  not  zero),  rather  than  maintaining  a  constant  value  across  the  beam 
as  is  assumed  theoretically. 

While  both  the  low  and  high  intensity  bright  screening  solitons  depend  on  the 
ratio  of  the  soliton  peak  intensity  to  the  background  plus  the  dark  intensity,  there  are  three 
striking  differences.  The  first  is  that  the  lowest  voltage  required  for  trapping  a  bright 
soliton  is  obtained  when  this  ratio  is  »  2.4  in  the  low  intensity  regime  and  5.5  in  the 
high  intensity  regime.  Indeed,  in  our  experiments  this  "most  favorable  point"  in  the  high 
intensity  regime  is  shifted  towards  higher  intensity  ratio  values  as  compared  to  the  low 
intensity  curve.  Second,  since  the  dark  irradiance  is  extremely  low  in  SBN,  i.e., 
milliwatts  per  square  centimeter  (or  less),  the  low  intensity  regime  can,  in  principle,  be 
reached  without  the  use  of  a  background  beam  (lb=0).  For  the  high  intensity  case, 
however,  the  background  beam  is  essential  and  is  many  orders  of  magnitude  higher  than 
the  dark  intensity.  Third,  the  slope  of  the  low  intensity  curve  is  significantly  greater  than 
the  corresponding  slope  for  the  high  intensity  regime  when  u^o  >  5.  This  difference  in 
slope,  for  example,  leads  to  a  soliton  width  in  the  high  intensity  regime  that  is 
substantially  lower  than  the  corresponding  soliton  width  in  the  low  intensity  regime  for  a 
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given  intensity  ratio.  As  seen  in  Fig.  2,  this  difference  in  slope  is  one  of  the  striking 
differentiating  features  supported  by  the  experimental  results  presented  in  this  paper. 
Since  precise  theoretical  predictions  only  exist  for  the  1-D  case,  a  1-D  experiment  was 
necessary  for  a  valid  comparison  with  theory.  However,  2-D  high  intensity  screening 
solitons  are  even  more  interesting,  both  for  fundamental  reasons  and  for  applications. 
Fig.  3  shows  typical  results  for  a  2-D  high  intensity  experiment  for  the  same  conditions  as 
the  1-D  experiment.  Although  not  shown  in  Fig.  2,  the  slope  of  the  experimental  2-D 
high  intensity  existence  curve,  when  u^o  >  5,  was  identical  to  the  1-D  high  intensity  case. 

In  conclusion,  we  have  demonstrated  that  screening  spatial  solitons  in  both  one 
and  two  dimensions  are  observable  using  8  nanosecond,  MW  optical  pulses.  The 
behavior  of  these  high  intensity  screening  solitons  differs  significantly  from  the 
previously  reported  low  intensity  c.w.  screening  spatial  solitons,  but  are  accurately 
described  by  their  existence  curve.  These  results,  will  now  be  combined  with  quasi-phase 
matching  to  generate  efficient  parametric  conversion  materials. 
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